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Abstract 

It is well-known that the number fluctuation in the grand canonical ensemble, 
which is directly proportional to the compressibility, diverges for an ideal bose 
gas as T ^ 0. We show that this divergence is removed when the atoms 
interact in one dimension through an inverse square two-body interaction. In 
two dimensions, similar results are obtained using a self-consistent Thomas- 
Fermi (TF) model for a repulsive zero-range interaction. Both models may 
be mapped on to a system of non-interacting particles obeying the Haldane- 
Wu exclusion statistics. We also calculate the number fluctuation from the 
ground state of the gas in these interacting models, and compare the grand 
canonical results with those obtained from the canonical ensemble. 
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I. INTRODUCTION 



Consider an ideal bose gas at low temperatures. In the grand canonical ensemble (GCE), 
the compressibility and the number fluctuation of this system diverge at low temperatures [|l|. 
To quote Landau and Lifshitz |]l|: 

...in a bose gas at temperature T < Tc, the pressure is independent of the volume, 
i.e. the compressibility becomes infinite. Accordingly... this would imply that 
the fluctuations of the number of particles also become infinite. This means 
that, in calculating fluctuations in a gas obeying bose statistics, the interactions 
between the particles cannot be neglected at low temperatures, however weak this 
interaction may be. When the interactions, which must exist in any actual gas, 
is taken into account, the resulting fluctuations are finite. 

One purpose of this paper is to demonstrate the validity of this statement in one and 
two dimensional traps for a special class of interactions. 

In one dimension, the quantum many-body problem of particles in a harmonic oscillator 
interacting with an inverse square two-body potential is exactly solvable @J^. Moreover, 
it is known that the global properties of these interacting bosons are the same as those 
of non-interacting particles obeying the Haldane-Wu generalized exclusion statistics (ideal 
haldons) [^,0. Using this mapping, we show that the number fluctuation is finite as T — > 0, 
no matter how weak the interaction strength is. In this model, since the exact correlation 
function is known for some specific strength parameters of the interaction it is also 
possible to verify explicitly the well-known relation between its integral and the number 
fluctuation. 

In two dimensions, there is no suitable exactly solvable model for our purpose. We 
therefore consider the mean-field model of bosons trapped in a harmonic oscillator, and 
interacting pair-wise with a zero-range repulsive pseudo-potential 0. In the absence of this 
interaction, the number fluctuation diverges at T = Tc. When the interaction is present, 
however weak, we show that the compressibility and hence the number fluctuation of the 
system is finite right down to T = 0. Moreover, they are shown to be identical to those of 
a collection of non-interacting haldons. 

Another quantity of some interest is the number fluctuation of particles from the ground 
state of the system, which is present even when the total number of particles in the trap is 
fixed. Consider a dilute gas of bosons in a trap at T = 0. The system is in its ground state. 
When a certain amount of excitation energy is given to the system, it may be absorbed 
in many possible ways, so that the number of bosons remaining in the ground state is not 
fixed. This number fluctuation for non-interacting bosons in a harmonic trap has been 
calculated by a number of authors as a function of the excitation energy or temperature. 
Thermodynamic identities and particle number fluctuations in weakly interacting BEC have 
also been analysed when the particle number is fixed For ideal bosons, the number 

fluctuation from the ground state diverges at low temperatures in GCE, but this can be 
avoided using more careful canonical, or microcanonical treatment. When the inverse square 
pair-wise interaction is used in one dimension, this divergence is again removed. Further, 
in this interacting model, we can also perform the canonical ensemble (CE) calculations, 
and compare with the GCE results. We find that even though the ground state number 
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fluctuation goes to zero as T — in both GCE and CE, the very low-temperature behaviors 
are rather different . At higher temperatures the results from CE and GCE tend to coincide. 
Similar studies can be made for fermionic systems at low temperatures. Since the interactions 
that we have used are of statistical character ||5|JTT[], our results interpolate from bosonic to 



fermionic behavior with the variation in the strength of the interaction. 

The outline of the paper is as follows: Sect. II deals with the total number fluctuation in 
GCE for interacting models. Sect. II (a) is devoted to a discussion of the fluctuations in one 
dimension, and II (b) to fluctuation in two dimensions. In Sect. Ill, the ground state number 
fluctuations for the one-dimensional interacting model in GCE and CE are calculated and 
compared. Particular attention is paid to the low temperature behavior (see also Appendix). 
We conclude the paper with a brief discussion of the results. 

II. FLUCTUATIONS IN GCE 

The number fluctuation in a gas in the GCE formalism is defined by {6N)'^ = (< A^^ > 
— < N >^), where the angular brackets denote ensemble averaging. For an ideal bose gas, 

oo 

(5A^)' = EK)(K) + i), (1) 

k=0 

where the single particle occupancy {uk) at a given temperature T and energy are given 
by the bose distribution function 

K) = rn. ^ ^1 T (2) 

exp[/3{ek - /i)J - 1 

Here f3 = 1/T with the Boltzmann constant = 1, and fi denotes the chemical potential. 
The isothermal compressibility Xt of ^ S^is of density po is related to the number fluctuation 
of the system in GCE : 

= TpoXt ■ (3) 

The problem with the number fluctuation in a bose gas in GCE is obvious from Eq. (p. 
If there is BEG, then a macroscopic fraction of the particles occupy the ground state for 
T < Tc, so that {SNqY ~ A^^. Even if there is no BEG, the same is the case as T — 0. In 
the thermodynamic limit, therefore, the fluctuation diverges below the critical temperature, 
or in any case at T = 0. This is manifestly incorrect, since at T = 0, all the bosons are in 
the ground state, and the number fluctuation should vanish. In the presence of interactions, 
however, this number fluctuation is expected to be finite. We demonstrate this by analyzing 
two models in one and two dimensions. 



A. Fluctuations in a one dimensional model 

We first consider the exactly solvable one dimensional Calogero-Sutherland Model (GSM) 
of a system of interacting particles with the Hamiltonian MM 
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N 



1=1 



^'52 1 2 2' 

H — muj X,- 



2m dxj 2 



,51 (4) 



with the dimensionless couphng parameter g > 0. The particles are confined in a harmonic 
well and the thermodynamic limit is obtained by taking u; — > as 00, with uN = 

constant. In the thermodynamic limit, the properties of the system are translationally 
invariant, and would be the same if the particles were on a line, or a circle, instead of a 
harmonic confinement. To make the problem well-defined quantum mechanically, we have 
to demand that the wave functions go to zero as |xj — Xj\^ whenever two particles i and j 
approach each other. Since the particles cannot cross each other, we may choose the wave 
function to be either symmetric (bosonic) or antisymmetric (fermionic). For g = and 1, 
the model describes free bosons and free fermions respectively. 

Using the relation between the integral of the correlation function and the number fiuctu- 
ation, we now show that the number fiuctuation vanishes at zero temperature in the above 
interacting model unlike the ideal boson result in GCE. If z/(r) denotes the two-particle 
ground state density- density correlation function in the ground state, with r = jxi — X2I, 
then the number fiuctuation is related to the correlation function as 

{6Nf 



N 

Note that the ground state correlation function z/(r) is defined only for r > 0. However, in 
computing the above integral it is necessary to assume z/(r) to be even function, and extend 



the domain of integration to negative values of r |12[. Unlike the one-particle off-diagonal 
density matrix, z/(r), by definition, is related to the diagonal element of the the two-particle 
density matrix, and is the same in GSM for bosons or fermions. In this section we work 
in the bosonic basis. The correlation functions are known exactly in the GSM for three 
values of g independent of whether the particles are bosons are fermions and are given (in 
the thermodynamic limit) by |^ 



= 1 : z/(r) = s(r) 



2 



sin(7rr)'' ^ 

vrr 
ds /■°° 



(6) 



(1 s 

g = l/2:u{r) = s{ry + -l dt [s{t)] (7) 



2 ds{2r) r2r 



g = 2: u{r) = s{2ry - dt [s{t)], (8) 

where the Fermi momentum kp is set equal to vr so that the maximum central density is 
unity. For all three forms of z/(r) given above, explicit calculations show that 

i-OO 

i^{x)dx = / [s{r)fdx = 1, (9) 



independent of the value of g. Substituting this result in Eq.(^ it follows that for interacting 
bosons in GSM the fiuctuation vanishes identically at T = 0. 

While we cannot obtain the exact z/(r) in GSM for all g, the same may be calculated for 
all values of g in the harmonic lattice approximation. The correlation function so obtained 
compares very well with the exact correlation functions for (7 = 1/2, 1, 2 and is given by ||13 
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where po = is the average density and 

, ^, ^, 1 /■'^ , 1 - cos (wn) 1 1 - cos (yn) ,,,, 

AFn,0 = — / dy ^ = — / dy 11 

ZTT^ Jo y — y /2tt ZTi^ Jo y 

The above expression is given for completeness and its exact form is not needed for further 
calculations. Again integrating over the real line we get a result identical to that obtained 
using the exact correlation functions in CSM. Thus the fluctuation vanishes identically for 
all g in this approximation at zero temperature. However, the result does not give any 
indication of the behavior of the fluctuation at finite temperature. To do this we take 
recourse to the mapping between the CSM and the exclusion statistics first proposed by 
Haldane through a generalized Pauli principle |5|,|17|. 



A crucial property of exclusion statistical interactions is that they should cause shifts in 
single particle energies at all scales |jT^ (see next section). This property is realized by a large 
class of one dimensional models of interacting fermions where Fermi liquid theory breaks 
down [1^,16]. In fact it has been shown exactly that quasiparticles with nontrivial exclusion 



statistics exist in a class of models that are solved by the Bethe ansatz |T^. Thus the result 
obtained below should be, in principle, valid in a large class of models with interactions. In 
particular, it is well known that the interacting particle of CSM may be regarded as ideal 
exclusion statistics particles or simply haldons 0. The thermodynamic properties of an ideal 
gas of exclusion statistics particles have been investigated widely 2 1 1 . The distribution 
function has been computed and is given by 

{w{e)+g) 

where w{e) is the solution of the equation 

w{e)3{l + w{e))'^^-^^ = e^(^-^) (13) 
At zero temperature we have, 

(n) = -, for efc < (14) 

g 

and zero otherwise. 

Note that the distribution function reduces to the usual Fermi and Bose distribution 
functions for g = 1 and g = respectively and in general g is regarded as the exclusion 
statistics parameter. Indeed, the statistical parameter g is precisely the interaction strength 
as given in Eq. (^ in one dimension. However, the distribution function as given above 
is valid in general and not necessarily restricted to one dimensional models. The following 
discussion is therefore used for illustration in the case of one-dimensional model but not 
restricted to this case alone. 

The number fluctuation at a given energy is given by. 
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(.«0^^r^. (15) 



Substituting for (n^) from Eq. (|I^), we have for total number fluctuation |TB 



(SNf = J2M'^+Wk){nkf (16) 

k=0 

oo 

= J2{nk){l-g{nk)){l + {l-g){nk)). 

k=0 

The number fluctuation vanishes at T ^ since — > 1/g below the Fermi energy and zero 
otherwise. This result holds no matter how weak the interaction strength is. However, at 
g = 0, the bosonic limit, the number fluctuation diverges as noted earlier. In this exactly 
solvable model, we have thus shown that interactions do remove the fluctuation catastrophe 
encountered in the ideal Bose gas. 

While these results have been derived in one dimension, extension to higher dimension 
is non-trivial since there is no suitable exactly solvable many body model. However, it has 
been shown that models with short range interactions in two dimension may be regarded 



as obeying exclusion statistics in the mean- field picture [Q. We discuss the fluctuation in 
these models in the next section. 



B. Fluctuations in a two dimensional model 

We consider a two-dimensional system of bosons interacting via a zero-range repulsive 
pseudo-potential. The quantum dynamics is then approximated by the following Hamilto- 
nian 

^4 (17) 

where the momenta and coordinates are planar vectors. The one-body potential generated 
by the above zero-range interaction (including exchange) is 

U{n{r)) = ^gn{r), g = 2~g , (18) 
m 

where ?T,(r) is the local number density of the system. In two- dimensions, g > plays the role 
of the statistical parameter, with (7 = for non-interacting bosons. At finite temperature, 
for T > Tc, the Thomas-Fermi approximation yields 



n r 



■• 2m 

where the Thomas-Fermi mean potential V{r) is given by 



V{r) = Voir) + g n{r) . (20) 

m 
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Here Vo(r) is the one-body harmonic trap. Note that Eq.(^9|) is vahd only in the absence of a 
condensate. It has been shown [^,0]. however, that for a nonzero positive g, a self-consistent 
solution of this equation satisfying / n{r)d?r = N may be obtained right down to T=0. This 
solution has a lower free energy than the one with a condensate [§], so we may take = 
for g > 0. 

The momentum integration may be done analytically: 

TTl 

n{r) = In [1 - exp[-/5(r(r) - /i)]] . (21) 

The local number fluctuation (between r and r + dr) in GCE is given by 

{SNf = T^, (22) 

where the coefficient of the temperature T on the rhs is related to the compressibility. Taking 
the derivative of the local density with respect to the chemical potential we have 

= exp[{V{r)-f,)P]-l + g ' ^^^^ 

Note that yU is a function of temperature, and is determined by the condition that 
/ n{r)(Pr = N, and in the thermodynamic limit it approaches the lowest energy eigen- 
state as the temperature goes to zero. A few remarks on the thermodynamic limit are in 
order: The thermodynamic limit is reached when N oo and uj ^ 0. However, in the 
limit of no confinement the density of states becomes a constant and there is no critical 
temperature below which condensation takes place. Preserving the density of states as in a 
two dimensional oscillator, the condensation temperature of an ideal bose gas is given by 

= (6/7r2)i/2 

Thus the limit — * oo and ^ is obtained keeping Tj°^ constant. In an ideal bose 
gas {g = 0), no self-consistent solution of Eq. (^) can be found for a fixed below this 
temperature Tj°) . However, when > no matter how small, the self-consistent solution of 



Eq. ( pID may then be found for all T > 0. 



In this limit the fluctuation is given by, 

(5iV)2 = , (24) 

2nh\^p[C-^gn-f^)P]-l + g 

where n is the constant density given in the thermodynamic limit. (For the one-dimensional 
case, this was denoted by po earlier.) 

In the absence of interaction, g = 0, the chemical potential /i goes to zero at T = 
and the above expression diverges as expected. However, as in the one dimensional case, 
the fluctuation remains flnite and approaches zero as T — when g is flnite, however weak 
the interaction may be. 

The same result may be obtained in the non-interacting exclusion statistics description 
in two- dimensions. The local density as a function of the radial coordinate is given by H], 
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n[r] = I , 25 

where the local variable w{p,r) is defined through Wu's equation within TF approximation 

2 

w^il + wy-'^ = exp[/?(|^ + Voir) - /.)], (26) 

and g is the exclusion statistics parameter which we will identify with the interaction strength 
in the mean-field picture. We have retained the trap potential Vq as in the interacting picture. 
Once again the momentum integration may be done easily and we obtain, 



n{r) = TT- In 



l+Wo 



(27) 



Wo 

where the local variable Wo{r) is determined through 

w'oil + Wo)'-' = explPiVoir) - /i)]. (28) 
The fiuctuation is then found by using Eq. (|22| ) and we have 

SN^ = tMI, (29) 

= T— 2 . (30) 

2TTh^ wo + g 

In the bosonic limit taking = in Eq. (|^), we have Wo = exp[/?(Vo — fi)] — l. Substituting 
this in the above expression for the fiuctuation in the thermodynamic limit it is easy to 
see that at Tc there is divergence. However, for positive definite g there is no divergence. 
Furthermore the equivalence between the non-interacting exclusion statistics picture and 
the mean-field description is established using the following relationship 

wo{r)=exp[P{V{r)-fi)]-l, (31) 

where V{r) is the self consistent mean-field potential. Substituting this in Eq.(pO|), and 
taking the thermodynamic limit yields Eq.(^4|). Note that this equivalence allows one to 
calculate fiuctuations in either the mean field picture or in the non-interacting exclusion 
statistics picture. Indeed this holds for the computation of other global thermodynamic 
quantities as well. 



III. GROUND STATE FLUCTUATION IN THE CANONICAL ENSEMBLE 

The total number fiuctuation defined in the previous section in the GCE is obtained by 
summing over all the single particle states. In the canonical ensemble however, the number 
of particles is fixed and therefore the fiuctuation in particle number has to be defined with 
respect to a reference state. One way defining the fiuctuation is to look at the ground state 
occupancy as a function of temperature, which is present even when the total number of 
particles in the trap is fixed. At T = all the particles are in the ground state. At a 
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nonzero temperature (or excitation energy), there are many ways of exciting the particles 
from the ground state, leading to a fluctuation in the ground state population. This number 
fluctuation for non-interacting bosons in a harmonic trap has been calculated by a number 
of authors as a function of the excitation energy or temperature. This has also been 
calculated for fermions in the CE by us in a previous publication Unlike the case of 
GCE for bosons, the ground state fluctuation in the CE is finite at all temperatures. 

We may extend the analysis of fluctuations in CE to particles interacting via the in- 
verse square pair-wise interaction in one dimension. It is more convenient to perform the 
calculation in the fermionic basis, although the formulae given here are applicable for both 
interacting fermions, or interacting bosons. In this section, quantities like energy and ground 
state number fluctuation of the interacting system with interaction strength g will be de- 
noted by a bracketed superscript, e.g. E^^\ (((5Aio)^)*-^^. We recall that the spectrum of the 
CSM Hamiltonian given in Eq. (P is exactly known. The states may be labeled by a set of 
fermionic occupation numbers {n^}, k = 1, oo, = 0, 1. The energy E^^^ of the system 
in the fermionic basis is then given by, 

E^'^ink} = E e^rik - c.(l - g) ^^^ ~ ^\ (32) 

k=l 

where = {k — ^hu) denotes the harmonic oscillator energy levels and = Z^fcLi'^fc- As 
can be seen from Eq. (|32D , the effect of the interaction is that each particle shifts the energy 
of every other particle by a constant huj{g — 1). The energy functional can also be written 
as 

oo oo 

E^^'^{nk} = 'Yeknk - uj{l - g) ^ ^k^^kz- (33) 

k=l ki<k2=l 

The exact spectrum of the model is thus reproduced by an effective Hamiltonian of quasi- 
particles with constant density of states and constant Landau parameters. As mentioned 
before, this scale invariant energy shift is the basic reason for the occurrence of nontrivial 
exclusion statistics where g plays the role of exclusion statistics parameter with = 0, 1 for 
bosons and fermions. 

The general canonical partition function in any basis is written in the occupation number 
representation as 

Zl^^ = J2 e-'3^'''^'^^>. (34) 

Using the energy spectrum in CSM given in Eq. (|3^) , the N-particle partition function in 
this one dimensional model is given by 

Zi^^ = e'^(-^)^Z-, (35) 

where (3 = jShuj and is the N particle fermion partition function. Setting g = 0, the 
bosonic partition function is obtained, 

= e^^Z-. (36) 
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Combining Eqs. (|35D and (|36|) we may write the partition function for CSM as, 



= {Z^YiZ^y-^. (37) 

The canonical partition function given above is exact in CE and may be used for calcu- 
lating the thermodynamic properties of the system in CSM within the canonical ensemble 
formalism. The moments of the occupation number are related to the partition function by 

^i-l; (-If)- 

where yk = exp(— /3efc). Therefore, it follows that 

{(fc.)^)<»' = y,?M^. (40) 

dyk 

Using Eqs. (p^)-(|3^), {nk)^g) can be expressed in terms of those of fermions and bosons via: 

{nk)^^^=g{nkf + {l-g){n,f (41) 

Unlike {nk)^^\ (^fc)''^'' does not have a simple form as Eq. (^Tj). However, the expression for 
the fluctuation in the occupation number does: 

{{5n,f)^g) = g{{5n,)Y + (1 - g)ii5n,)y. (42) 

Eq. (|42|) gives only the fluctuation in the occupation of a given level k, while the quantity 
we are seeking is the ground state number fluctuation. The latter is formally deflned in any 
ensemble as: 

= Y.^6n,f = Y.{{ni) - {n,f) (43) 

k k 

where the sum k runs over only the levels which are completely occupied at zero temperature. 
Thus, in an ab initio calculation, one would formally sum over the quasiparticle levels which 
are occupied at T = to get {{SNof'^g). Fig. 1 shows the level flow in CSM as a function 
of g obtained from Eq. (|33D at T = 0. It can be seen that as g changes from the fermionic 
to the bosonic end, the number of levels contributing to the ground state remains constant, 
while the Fermi energy decreases accordingly. This means that one may obtain ((5A^o)^)^^^ 
by simply substituting the ground state fluctuations for fermions and bosons, ie: 

{{SNofY^^ = g{{5No)Y + (1 - g)mo)Y (44) 
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IV. DISCUSSIONS 



In the previous section we have discussed the ground state number fluctuation in GCE 
and CE in an interacting system. We now compare the results obtained in these two en- 
sembles. Due to numerical difficulty at low temperatures, we are able to do the comparison 
only for = 10 particles. Despite this limitation, some interesting points can be made. In 
Fig. 2(a), we show the behavior of the relative ground state fluctuation against temperature 
of interacting bosons for both ensembles for interacting strengths of g = 0, 1/2, 1. The g = 
case corresponds to free bosons. Fig. 2(b) shows the low temperature region of the curves 
in Fig. 2(a). Note that as T ^ the GCE fluctuation for free bosons diverges as expected, 
whereas those for interacting bosons remain finite and approach zero. Note also that in the 
ideal Haldane gas picture the g = 1/2 case corresponds to semions and the g = 1 case is the 
non-interacting fermionic limit. The GCE fluctuation for semions is found using Eqs. (0), 
(P!B|) and (|T7P) where is determined by the constraint that the average total number of 
particles is N. 

Unlike in GCE, the CE ground state fluctuation for free bosons remain finite as T ^ 0. 
The ground state fluctuation of free bosons, however, approaches zero exponentially (see 
Fig. 2(b)), contrary to previous results that find a linear dependence with T all the way 
to T = [P|,P3|. In Appendix A we give the low temperature expansions of the fluctuation 
squared (SNq)'^ in powers of x, where x = . Clearly, at very low temperatures (^A'o)^ is 
independent of the number of particles A^ (in both GCE and CE). Therefore, the exponential 
behaviour of the fluctuation of free bosons at low temperatures should remain valid even 
in the large A^ limit. In the case of fermions the CE and GCE fluctuations are similar 
except at very low temperatures. Clearly, the CE and GCE curves of fermions approach 
zero differently as T ^ 0. For high temperatures, on the other hand, GCE and CE give 
identical results, as expected. 

It is to be noted from Fig. 2 (a) that the number fluctuation for free bosons and free 
fermions cross at a certain temperature, with the fermion one getting bigger at higher 
temperatures. This is to be expected, since the number of possibilities of creating holes 
within the fermi sea, and distributing particles above, increases more rapidly than for bosons, 
whose ground state has only one level. Semions {g = 1/2), which may be regarded either as 
interacting bosons, or as interacting fermions, have fluctuations intermediate between those 
of bosons {g = 0) and fermions {g = 1). 

In one dimension, we have calculated the number fluctuation of bosons interacting via 
GSM using the ideal Haldane picture for both GCE and CE. The next challenge is to find the 
exact fluctuation using combinatorics. This has been done for free fermions which correspond 
to interacting bosons with interacting strength of g = 1 p3|]. Work is in progress to calculate 
this exact fluctuation for a general g. 
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APPENDIX: LOW TEMPERATURE EXPANSIONS IN GCE AND CE 



The low temperature behaviour of thermodynamic quantities in GCE are well known 
for bosonic systems. However, a comparison between GCE and CE calculations at low 
temperatures is not usually discussed for either Bose or Fermi systems. Further, making use 
of some asymptotic expansions, the CE fluctuation for bosons was earlier found to be linear 
right down to T = 0. However, we give here the expansion of the fluctuation squared at low 
temperature in power of x, where x = e^, and show that the CE fluctuation of bosons is in 
fact exponential at very low T. In GCE only expansion for fermions is possible, since the 
fluctuation tends to infinity at low temperature for bosons. Both expansions are possible in 
CE. 

a. Grand Canonical Ensemble: In GCE the (fermionic) occupation number is: 

^^'^^^^^ = a:(>^-l) + 1 = a:^,~k-\/2) + i (45) 
for a one dimensional system. The ground state number fluctuation squared is given by: 

where kp is the Fermi level. At low temperatures, for the one- dimensional harmonic oscil- 
lator, jj, ^ jj,Q = N . Therefore, 

{{5NQf)GCE = ^/x-2x + Ax^''^ -Ax^ + Qx^'^ -Sx^ + %x'"^ -Sx\.. , iV > 4. (47) 

Note that the first few terms are independent of A^. 

h. Canonical Ensemble: In CE the first and second moments of the occupation number 
are known [24|: 

1 ^ 

K) = ^E(±)'^V^'=^^-.> (48) 

1 ^ 

= ^ E(±)''^' ± (j - 1)] (49) 

where the upper and lower signs refer to bosons and fermions respectively. Summing over 
the ground states up the Fermi level gives the fermionic ground state number: 

(A^o)cE = ^ Y.{-^r'^'''^-^ Zn-,. (50) 

Zjv J V — X-i 

where we let = /c — 1/2, with k = 1,2, .... Therefore, 

AT N 
k=l k=l 

= x + 2x^ + ... , N>4:, (51) 
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where we have used (uk) — (rip.) and again the first few terms are independent of the system 
size A^. 

For bosons the ground state consists of one single lowest level, the low temperature 
expansion of the number fluctuation is given by 

{{SNor)cE = {nl) - (no)' 

= x + 3x2 + 4x^ + 7a;^ + ... , N > 4. (52) 

Again as in the fermionic case the first few terms in the low temperature expansions are 
independent of the system size. Indeed it is interesting to note that in CE, the fluctuations 
in both the systems approach zero as T — > in exactly identical fashion. 
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FIGURES 



FIG. 1. The level flow of quasi-particle energy levels in CSM as a function of g. 

FIG. 2. (a) The ground state fluctuation for the one-dimensional CSM system in GCE {Eq. ( [T7| ) 
and CE {Eq. ( ^4|) as a function of temperature for N = 10. We show the results for fermions, 
bosons and also semions {g = 1/2). 

FIG. 3. (b) Same as in (a) but using the low temperature expansions given in the Appendix 
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